Crystalline phases in chiral ferromagnets: Destabilization of helical order 
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In chiral ferromagnets, weak spin-orbit interactions twists the ferromagnetic order into spirals 
leading to helical order. We investigate an extended Ginzburg-Landau theory of such systems 
where the helical order is destabilized in favour of crystalline phases. These crystalline phases are 
based on periodic arrangements of double-twist cylinders and are strongly reminiscent of blue phases 
in liquid crystals. We discuss the relevance of such blue phases for the phase diagram of the chiral 
ferromagnet MnSi. 
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I. INTRODUCTION 

In metallic magnets without inversion symmetry - so- 
called chiral ferromagnets - spin-orbit coupUng effects de- 
termine the physics decisively: they can twist the magne- 
tization into a helix in the ordered phase of the material. 
Prominent examples are MnSi and FeGe, which have for 
a long time been known to exhibit spiral orderi. Re- 
cently, interest in such materials has been renewed with 
a number of experiments that on the one hand found 
non-Fermi liquid behaviour in a large temperature and 
pressure region of the phase diagram^, and on the other 
hand uncovered signs of a peculiar partially ordered state 
in neutron scattering experiments"^. 

In the ordered state of MnSi, helical order is observed 
in neutron scattering experiments in the form of Bragg 
peaks situated on a sphere in reciprocal space. The ra- 
dius of this sphere is proportional to the inverse pitch 
of the helix, and the peaks are positioned in the (111)- 
directions. This locking of the helices is due to higher 
order spin-orbit coupling effects. For temperatures down 
to 12 K the phase diagram exhibits a second-order phase 
transition out of the helical into the disordered, isotropic 
phase when external pressure is applied. Below 12 K, 
helical order is lost at external pressure of 12-14 kbar via 
a first-order phase transition into a "partially-ordered" 
state. This state seems to retain remnants of helical 
order on intermediate time and length scales: neutron 
scattering experiments reveal a signal on the surface of 
a sphere in reciprocal space, which is resolution-limited 
in the radial direction but smeared out over the surface 
of the sphere^. Maxima of the signal now point into the 
(llO)-directions but are no longer sharply peaked. 

Motivated by these experiments, we investigated ways 
of destabilizing helical order in chiral ferromagnets and 
the possible new phases that can exist besides helical 
order. Such phases are well known for cholesteric liquid 
crystals^, i.e. chiral nematics, where several so-called blue 
phases have been observed. In these phases (more pre- 
cisely in the phases I and II) crystalline order is formed 
which can be interpreted as a periodic network of ordered 



cylinders, see below. The lattice spacing is often of the 
order of a few hundred nanometers, leading to the color- 
ful appearance of these phases (including the color blue 
in some variants). 

As the Ginzburg-Landau theory for the director of a 
chiral nematic liquid is very similar to that for the vec- 
tor order parameter of a chiral magnet (see discussion 
below), the question arises whether similar phases are 
realized in magnetic systems. This question has been 
addressed in some detail by Wright and Mermini many 
years ago. They showed that amplitude fluctuations, 
which are essential for the stability of blue phases, cost a 
factor 3 more energy for ferromagnets compared to liquid 
crystals and concluded that such phases do not appear 
within a Ginzburg Landau theory with local interactions. 

Motivated by the physics of MnSi, several groups have 
reexamined this issue. RoCler et ali^ added a further pa- 
rameter to the Ginzburg Landau theory, allowing them 
to reduce the energy of amplitude fluctuations arguing 
that such a term might arise from higher-order fluctu- 
ation corrections. The term considered in Ref. Q is, 
however, non-analytic in the Ginzburg-Landau order pa- 
rameter. In the presence of such a term, they were able to 
show that a crystalline array of cylinders can have lower 
energy than a uniform helix. We will discuss a similar 
structure below. An alternative route (followed also by 
us) to destabilize the helical solution are non-local inter- 
actions as suggested by Binz et oIf . They considered 
crystals formed from superpositions of helices to be re- 
sponsible for the partially ordered phase of MnSi. 

As pointed out by Wright and Mermini, is is useful to 
distinguish two rather different limits when investigating 
the physics of blue phases. Here one has to compare 
the correlation length ^ for amplitude fluctuations of the 
ordered state (i.e. the width of a typical domain wall) 
with the wavelength 27r/go of the helical state. For ^go ^ 
1, the high-chirality limit, the magnetic structure can 
be described by a superposition of a few helices. This 
approach was taken by Binz et al^. 

The factor ^qq is proportional to the strength of spin- 
orbit coupling! and therefore expected to become large 
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only very close to the second-order phase transition 
where ^ diverges. As the pitch of the helix'^ in MnSi, 
is very large compared to all other mi- 
croscopic length scales, we think that the low-chirality 
limit, ^go ^ 1, is more appropriate for the description 
of the high-pressure phase of this material. In this limit, 
amplitude fluctuations of the order parameter cost more 
energy than twists of the phase and therefore one has 
to look for order-parameter configurations with as little 
amplitude fluctuations as possible. 

In the following we will first investigate which terms in 
the Ginzburg Landau theory destabilize the helical state. 
Then we will suggest variational solutions appropriate in 
the low-chirality limit. Finally we investigate experimen- 
tal consequences of the resulting structures. 

II. CHIRAL MAGNETS AND BLUE PHASES 

Starting point of our investigation is a Ginzburg- 
Landau theory at finite temperature, assuming that all 
modes with non-zero Matsubara frequencies are massive 
and can be integrated out. 

Up to second order in spin-orbit coupling, the 
Ginzburg-Landau theory for chiral ferromagnets is given 

/(^) = f EC^^-^^)' + 7M • (V X M) + /fm, (1) 

where M = M{r) is the position-dependent magneti- 
zation and /fm = |(M)^ -f- u{M)^ is the Landau free 
energy of the underlying ferromagnet. Spin-orbit cou- 
pling is present in the form of the Dzyaloshinsky-Moriya 
interaction oc 7. In order to facilitate the following dis- 
cussion, we also provide a second form for ([T]). It can be 
rewritten by setting M — Xn, where A is the amplitude 
and n is the direction of the magnetization M: 

fir) = f (V.n, + ls.,,n,y + ^(VA)^ - ^A^ 

+ /fm(A). (2) 

with /fm(A) = |A2 -hwA". 

Terms that break rotational symmetry are of higher 
order in spin-orbit coupling and take the form e.g. 

Bi [{dM^ + [dyMyf + (9,Af,)2] 

+ B2 [{dlMf + {dlMf + {dlMf] 

+ B,{M^ + M^ + M^). (3) 

(to fourth order in spin-orbit coupling) for the point 
group P2i3 of MnS:^. These terms will be neglected for 
the moment. 
The helical state 

^hohx^^-j ^ xcos{qaz) -f y sin(goz), (4) 

with go = 7/0^ can be shown to be the lowest energy 
state of ([IJ if one assumes that the amplitude A of the 



magnetization is constant and was argued in Ref. [4| to be 
the only ordered phase possible for chiral ferromagnets. 
Higher order terms in the Ginzburg-Landau expansion 
can, however, serve to destabilize this helical order and 
induce other phases. The simplest (i.e. of lowest order 
in spin-orbit coupling) of these terms will be the focus of 
our investigation in this paper: 

(^M,VM}j =liVMy. (5) 

This term acts only on the amplitude of the magnetiza- 
tion, not on its direction. Therefore, it gives no contribu- 
tion to the free energy density in the helical phase, which 
has a uniform amplitude. If ^ has a negative prefactor, 
then in a certain parameter regime, it can be expected to 
destabilize helical order in favour of an order parameter 
with a fluctuating amplitude. In the following we will 
use the shorthand (Af VAf )^ for the expression ([5]). 

In order to determine the new states that could be 
stabilized by (O, it is instructive to explore the close 
analogies between chiral ferromagnets and chiral liquid 
crystals^. The order parameter of chiral liquid crystals is 
a director, and as a consequence topological defects are 
fundamentally different in both systems. However, chiral 
liquid crystals can be described by a free energy density 
that is quite similar in form to the only differences 
being an additional cubic term in A and an extra factor 
1/3 in front of the term cx (VA)^. 

Locally, the blue phases are based on a configuration 
of the magnetization that can be shown to be even lower 
in energy than the helix: 

n^*(r) = i:cos(gr) — 0sin(gr), (6) 

in cylinder coordinates, see Fig. [1] This can easily be 
seen4 by comparing the energy density of the uniform 

helix, — ^A^-f /fm(A), to the result obtained by plugging 
Eq. ^ into Eq. For A = const, and q — qq, the 

first two terms in ^ vanish for r = 0, and therefore the 
energy density at r = is given by — ^A^ -I- /fm(A), i.e. 

lowered by an extra factor — ^A^. 

This so-called "double-twist" configuration is cylindri- 
cally symmetric: sheets of constant magnetization are 
rolled up around a common cylinder axis (see Fig. [1]). 
The configuration ^ is, however, only favoured in the 
vicinity of the cylinder axis - the energy difference be- 
tween the helix and the double-twist configuration dimin- 
ishes and even becomes positive as the distance from the 
cylinder axis is increased. Isolated double-twist configu- 
rations therefore necessarily have to occur in the form of 
cylinders with an amplitude that becomes zero at a cer- 
tain distance from the cylinder axis. However, amplitude 
fluctuations cost energy, as can be seen from the second 
term in ([2]). In the case of liquid crystals, the prefactor 
of (VA)^ is small enough for double-twist cylinders to 
become energetically lower in energy than the helix (see 
discussion above). Crystals made of these double twist 
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netization are wrapped as cylinders around a common axis. 1 2 3 4 j- 



cylinders are indeed believed to lie at the heart of the 
blue phases I and II in liquid crystals^. In the case of the 
ferroniagnet, however, the energy cost of amplitude fluc- 
tuations outweighs the energy gains due to directional 
fluctuations within the double-twist structure, and for 
a free energy density given by ^ (single) double-twist 
cylinders do not occur-^. 

Adding the term ^ to the free energy density can 
invalidate this conclusion: if the prefactor of ^ is al- 
lowed to become negative, it can reduce the cost of am- 
plitude fluctuations and allow for the appearance of blue 
phases even in chiral ferromagnets. If ([5]) has a nega- 
tive prefactor, then the inclusion of higher order terms in 
the Ginzburg-Landau theory becomes necessary in order 
to obtain a stable solution and avoid e.g. an unbounded 
magnetization. 

In the following, we want to analyze the possible occur- 
rence of the analogue of blue phases in chiral ferromag- 
nets. By rescaling the magnetization, the free energy 
density and the momenta, the free energy density can be 
cast into the following form: 

/(r) = 5M^ + ^{VA'Uf + M • (V X M) + 

+ ^{MVMf +r]h,iM), (7) 

where C < Oi rj > and hi{M) is a term containing k 
powers of M and I derivatives (fc > 4, ^ > 2), which is 
used to stabilize solutions against an unbounded mag- 
netization and oscillations thereof. Possible choices are 
for example hi{M) = (M VM)^ ^(VMi)^, /i2(M) = 
(MVM)4 and hjiM) = M^{MVMf. In the rest of 
this paper, we will use hi{M) exclusively: this term is 
the one with the least number of derivatives and powers 
of M that on the one hand stabilizes single double twist 
cylinders and on the other hand is identically zero for the 
single helix. 

III. FROM CYLINDERS TO CRYSTALS 

As a first step, we investigated single double twist 
cylinders for S = 1/4. At this point, the helical and 
isotropic phase are degenerate with jhchcai _ j^isotropic _ 
0. If a single cylinder can now be shown to have negative 



FIG. 2: Amplitude function A(r) which minimizes ((Tj using 
the directional dependence (0), the stabilizing term hi (with 
^ — —6.5, -q — 5) and the boundary condition A(r) = for 
r > R. Here q — O.SSqo and R = 3.98 are free variational 
parameters where qo = 1/2 is the helix wave vector. The 
resulting free energy gain per length is given by J'^^' — —0.21. 



free energy, then the system can lower its free energy fur- 
ther by creating extended networks of double-twist cylin- 
ders: crystalline phases, the analogy to blue phases in 
liquid crystals, can be expected to form in chiral ferro- 
magnets. 

Setting M ri'"(r)A(r), with n'^'(r) given by ©, we 
calculated numerically the amplitude function A(r) that 
minimizes the free energy density ([7]) (for certain values 
of ^ and 77), subject to the condition that the magneti- 
zation drops to zero at a certain distance from the cylin- 
der axis. For a given 77 and sufficiently negative values 
of ^, single cylinders are stable configurations within a 
Ginzburg-Landau theory of the form ([7]) , see Fig. [2l In 
order to minimize its free energy, the system will try to 
produce many such cylinders, packed as tightly as possi- 
ble. In these configurations, the magnetization only has 
to drop to zero on lines or points, if at all. Considering 
that it is the competition between phase and amplitude 
fiuctuations that either stabilizes double-twist structures 
or not, it can be expected that crystalline arrangements 
of the helices can exist even in parameter regimes where 
single cylinders are unstable. 

The possible crystalline structures are subject to the 
condition that the magnetization matches where cylin- 
ders touch, in order to avoid discontinuities in the magne- 
tization. This condition is much more restrictive for fer- 
romagnets compared to liquid crystals^. We have found^ 
only two allowed structures constructed as networks of 
double-twist cylinders: a square and a cubic lattice, see 
Fig. [3) The square lattice is invariant with respect to ro- 
tations by multiples of 7r/2 around the z-axis, translations 
along the z-axis as well as a translation by a(i, ^,0) com- 
bined with time reversal {M — > —M). In this structure, 
the magnetization has to go to zero on the line (|-, 0, z) 
and symmetry-equivalent lines. The symmetry transfor- 
mations that leave the second structure (cubic lattice) 
invariant are cyclic permutations of the axes and a trans- 
lation by 0(5, 5, \ ) combined with time reversal. In this 
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FIG. 3: Crystalline structures built from double- twist cylin- 
ders: (a) Unit cell of the square lattice of double-twist cylin- 
ders in the x-y-plane. The magnetization on the cylinder axes 
as well as at the points where cylinders touch is shown explic- 
itly, (b) Unit cell of the cubic lattice of double-twist cylinders. 
At the points where cylinders touch, the magnetization has 
twisted to an angle of 45° from the cylinder axis. 
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FIG. 4: Phase diagram for r] = 0.05. In addition to the 
helical and the isotropic phases there is a parameter range 
where the cubic lattice of double-twist cylinders minimizes 
the free energy. 



structure, the magnetization vanishes only at a (f i f i f )i 
a(|,|,|), and symmetry-equivalent points in the unit 
cell. 

We calculated the free energy density of these struc- 
tures by means of a variational ansatz for the amplitude 
of a single double- twist cylinder: 

A(r)=yo-(ro-r)e-^/''^e(ro-r), (8) 

where j/oi '"i the double- twist wave vector q are vari- 
ational parameters. This ansatz is based on the numerics 
for a single cylinder. The cylinders were then arranged 
in crystals as shown in Fig. [31 The cutoff tq « -s/Stt/ {2q) 
was chosen for convenience to avoid within the calcula- 
tion the tiny overlap of cylinders which are far apart. 

A phase diagram as a function of the remaining three 
free parameters can now be computed. Here, we set rj 
to a fixed value, assuming that it is the least susceptible 
to variations of external pressure, and established the 
phase diagram as a function of the two remaining free 
parameters. The result for 77 — 0.05 can be seen in Fig.m 

One immediately notices that there is no parameter 
regime in which the square lattice shown in Fig. [3] (a) 
has the lowest free energy density; this is also true for 
different values of 77 not shown here. While the filling 
fraction of double-twist regions is certainly larger in the 
square lattice than in the cubic structure, the magnetiza- 
tion in the cubic crystal only has to go to zero at points 
and not on lines, as it is the case for the square lat- 
tice. Furthermore, in the cubic lattice the magnetization 
only has to twist outwards by 45 degrees until cylinders 
touch, compared to 90 degrees in the case of the square 
lattice (see Fig. ^ . These factors conspire to make the 
free energy of the cubic structure even lower than that of 
the square lattice. Within our variational ansatz we find 
that the transition from the crystalline to the isotropic 
phase is very weakly first order [the relevant free energy 
differences are less than 0.02 in units of Eq. ([7])]. Note, 
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FIG. 5: Ratio of double twist wave vector q [see Eq. ((Gjl] 
versus helix wave vector go [see Eq. Q] for rj = 0.05: q/qo 
varies considerably but remains of order 1 in the parameter 
region where the cubic crystal is stable. 

however, that our analysis might break down close to the 
transition, as one enters the high-chirality regime, see 
above. 

Since the double twist wave vector q is now also a vari- 
ational parameter, it is no longer necessarily identical to 
the helix wave vector go- In Fig- [H we show that q de- 
pends considerably on the microscopic parameters but 
remains of order qo . 

What signal can these structures be expected to pro- 
duce in neutron scattering experiments? As for any crys- 
tal, peaks in neutron scattering originate from the lat- 
tice structure itself, and the configuration of the magne- 
tization within the Wigner-Seitz cell only enters in the 
shape of form factors. The lattice constants are func- 
tions of the double-twist wave vector and determined by 
the requirement that the magnetization has to match 
where cylinders touch. For the square lattice, one ob- 
tains a = ^pl-n I q^ and for the cubic lattice a = 27r/q. 
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In neutron scattering, higher order Bragg peaks can 
also be expected to be be generated by the lattice struc- 
tures. The elastic scattering crossection^ can be calcu- 
lated from a Fourier transform of the magnetization: 

(^^) ^<^\kx{MiK)xk)\'. (9) 

In order to compare with experimental data, the quanti- 
ties of interest are the positions of the Bragg peaks and 
the relative intensities of the Bragg peaks. 

For the square lattice, normalizing intensities to give 
unity for the first reflection with Miller indices (10), the 
(21) and the (30)-reflections have intensities 0.08 and 
0.01, respectively. The invariance of the square lat- 
tice with respect to a translation by a(^, ^) combined 
with M —M constrains all Bragg peaks {hk) with 
h + k — 2n to vanish. In the case of the square lattice, 
the lowest order reflexes already account for 84% of the 
total intensity. 

For the cubic lattice, assuming that the lowest Bragg 
peak, i.e. the (lOO)-peak, has intensity 1.0, the (210)- 
and (300)-peaks have relative intensities 0.17 and 0.03. 
The (110)-, and (200)-peaks vanish for symmetry rea- 
sons, while the (lll)-peak vanishes as a consequence of 
our ansatz based on a linear combination of cylinders in 
X, y and z direction, and should really be non-zero as 
allowed by symmetry. Our approach suggests, however, 
that these peaks have small weight. The (lOO)-peaks of 
the cubic lattice only represent 46% of the total scatter- 
ing intensity /cubic- The (100)-, (210)- and (300)-peaks 
of the cubic lattice add up to 79% of /cubic- 

Anisotropic terms that orient the helix also act on 
the crystalline structures and determine their orientation 
with respect to the crystal lattice of the substance. For 
the square lattice, where all cylinders are arranged in par- 
allel, we find using Eq. ^ for Bi = B2 = 0, B3 ^ that 
weak anisotropic terms align the cylinder axes parallel 
to the preferred direction of the helix vector, i.e. either 
in the (111) or (100) direction, depending on the sign of 
B^. However, the orientation of the crystal in perpen- 
dicular direction is not affected to leading order by the 
anisotropy term, but higher-order terms would lock a per- 
fect crystal. Motivated by the "partial order" observed 
in the high-pressure phase of MnSi (see introduction), we 
investigate the expected signature in neutron scattering 
assuming that these higher-order terms are not effective. 
In this case, the square lattice will produce rings in planes 
normal to the (lll)-direction (the orientation of the helix 
in the low-pressure phase). These rings intersect to pro- 
duce maxima in the (llO)-direction on a circle with radius 
^/2q in reciprocal space. For parameters with ^/2q w qo, 
this is consistent with the observed signatures'^ in neu- 
tron scattering. Note, however, that at least within our 
model, the square lattice never has the lowest energy. 

In Ref. , it was argued that an amorphous texture 
of parallel cylinders (which the authors called skyrmions) 
aligned preferentially along the (111) direction would also 
produces such rings. However, such a scenario would not 



explain the resolution limited width in radial direction, 
i.e. at least on length scales of 2000A the square crystal 
would have to remain intact. 

In the case of the cubic lattice, the numerical calcu- 
lation of the free energy density shows that for B3 > 
the anisotropic terms are minimized if one of the axes of 
the cubic lattice is IKl/v^, 1/V2, 0) and the other axes 
are ||(-l/2, 1/2, 1/V5) and |j(l/2, -1/2, 1/^/2), respec- 
tively. While such structures would produce peaks in the 
(110) direction (as observed in the high-pressure phase 
of MnSi), one expects also considerable intensity rather 
close to (111) (as e.g. (-1/2,1/2,1/^2) differs only 
by about 10 degrees from (—111)). Experimentally, how- 
ever, one observes a minimum of the intensity in the (111) 
direction. For Bx = B2 = 0, B3 < 0, the system would 
prefer to align the cubic structure with the crystal lattice 
of MnSi. 

Binz et al.— have argued that the Bragg peaks in neu- 
tron scattering should be smeared out because the mag- 
netization has a varying amplitude and is therefore sus- 
ceptible to interactions with non-magnetic impurities. 
The same argument would apply for the crystalline struc- 
tures presented in this paper. 



IV. CONCLUSIONS 

In conclusion, we have constructed "blue phases" in 
chiral ferromagnets. Is it likely that these blue phases 
are realized in MnSi? The signatures in neutron scat- 
tering seem to be consistent with the square structure 
(also considered in Ref. 0), which is, however, never 
the ground state within the models which we considered. 
Our results suggest that the intensity in neutron scatter- 
ing is located on rings intersecting in the (llO)-direction. 
It would be interesting to check experimentally whether 
the intensity distribution arising from such a picture fits 
quantitatively. 

In the case of the cubic lattice, the positions of the 
expected maxima do not match the ones observed. In 
any case, the smoking gun experiment to detect crys- 
talline structures, would be the observation of higher or- 
der Bragg peaks (this applies also for other structures, 
e.g. those suggested in Refs. [1@|) and we hope that our 
estimates may guide future experiments in this direction. 

However, it is still unclear whether the partially or- 
dered state of MnSi is indeed a separate phase. Re- 
cent measurements of the thermal expansion coefficient 
give no trace of a phase transition from the partially or- 
dered state to the isotropic phase^^. Furthermore, /iSR 
experimentsii give no evidence for static order in this 
regime and it remains unclear on which time scales the 
partial order survives'^. 

It therefore remains an open question whether distinct 
phases other than the helical one are present in chiral fer- 
romagnets that are currently investigated experimentally. 
In fact, the partially ordered phase in MnSi might be 
more analogous to the blue phase III than to blue phases 
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I and II. While the precise structure of the blue phase 
III is not completely clear, it seems to show some type of 
crystalline short range order but is a liquid on long length 
scalesi^i^. Therefore it has been conjectured^'' that an 
amorphous arrangement of double-twist cylinders is rele- 
vant for this blue phase III in liquid crystals - this might 
also be true^i^ii^ for the partially ordered phase in MnSi. 
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